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We report the binomial multiplicative model for low impact energy fragmentation. 
Impact fragmentation experiments were performed for low impact energy region, 
and it was found that the weighted mean mass is scaled by the pseudo control 
parameter multiplicity. We revealed that the power of this scaling is a non-integer 
(fractal) value and has a multi-scaling property. This multi-scaling can be inter- 
preted by a binomial multiplicative (simple biased cascade) model. Although the 
model cannot explain the power-law of fragment-mass cumulative distribution in 
fully fragmented states, it can produce the multi-scaling exponents that agree with 
experimental results well. Other models for fragmentation phenomena were also 
analyzed and compared with our model. 
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1 Introduction 

Impact fragmentation of brittle solids is a typical nonlinear phenomenon. Small 
impact cannot make brittle solids cleave. However, large impact produces cracks 
irreversibly and makes brittle solids fissure to small pieces of fragments. This ubiq- 
uitous phenomenon can be seen even in our everyday lives. Thus, many scientists 
and engineers have studied this issue. As known well, cumulative distribution of 
fragment mass shows power- lawffl Oddershede et al0and Meibom and Balslev^l in- 
vestigated what controls the exponent of power-law distribution. They found that 
the exponent is determined by the dimensionality of fractured object. Ishii and 
Matsushita performed the 1-dimensional fragmentation experiments with long thin 
glass rodfl They dropped the glass rods from various heights. The cumulative dis- 
tribution obeyed power-law form at high dropping height, and obeyed log-normal 
form at low one. 

Recently, Kun and Herrmann investigated the damage-fragmentation transition 
for low impact energy collision by numerical simulation^. They used the granular 
solid disks colliding by a point to each otheJSI. The transition from damaged state to 
fragmentation state was observed by increasing the relative collision speed. They 
measured the critical exponents of damage-fragmentation transition and realized 
that scaling-laws of the percolation universality are satisfied near this transition 
region. On the other hand, Astrom et al. studied the low energy fragmentation 
using the random distorted lattice with elastic beam model and fluid MD model 
with LJ pair potential. They corrected that the critical behavior for low energy 
fragmentation differs from that of percolation. Oddershede et al. said the fragmen- 
tation process is a kind of self-organized critical phenomenon^. However, most of 



BinomialModel: submitted to World Scientific on February 2, 2008 1 



experiments examined only on high imparted energy fragmentation. There are no 
experiments on critical behavior of fragmentation by low inparted energy. 

In order to study the critical fragmentation, we performed simple experiments of 
fragmentation. We considered a simple binomial multiplicative scenario of critical 
fragmentation. In this paper, we report on results of detailed analyses on the model. 
In the next section, experimental results are presented. In Sec. |3 we introduce a 
binomial multiplicative model and analyze it. In Sec. 01 the results are compared 
with other possible models. 

2 Experiment 

We used glass tube samples as 2-D fractured objects. The tube was put between a 
stainless stage and a stainless plate. A brass weight was dropped to the stainless 
plate. The falling height was controlled on slightly higher than the point at which 
samples did not fracture. After fragmentation, we collected fragments and measured 
the mass of each fragment with an electronic balance. Fractured tubes have the 
approximate 2-D geometry (50 mm outside diameter, 2 mm thick, and 50, 100, 150 
mm length) . More detailed experimental setups are described in Ref. |H1 

According to Kun and Herrmann's result, the control parameter should be the 
imparted energy per unit sample mass e, and the order parameter should be the 
maximum fragment mass M max ^. The e is calculated as e = M w ghjM \ >1 where 
M w , g, h, and M b correspond to the mass of falling weight, the gravitational ac- 
celeration, the height of falling weight, and the mass of target sample, respectively. 
The log-log plot of maximum fragment mass M max vs. imparted energy per unit 
sample mass e is shown in Fig. ^ As can be seen in Fig. ^ Mnax and e relate with 
negative correlation, qualitatively. However, since the data in Fig.^contain large 
uncertainties, we cannot discuss quantitatively on the critical scaling by this plot. 
Therefore, we have to use another parameter to analyze quantitatively. 

Campi proposed a pseudo control parameter multiplicity /i in Ref. The \i is 
defined as, 



Where m m ; n , Mo, and M\ correspond to the smallest limit of fragment mass (we fix 
it at 0.01 g), the total number of fragments, and the total mass of the all fragments, 
respectively. The fragmentation critical point corresponds to the value fi = by 
this definition. In general, we can introduce the k-th order moment of fragment 
mass distribution M& as, 



where n(m) is the number of fragments of mass m. Certainly, Mq and Mi in Eq. 
JTJ are specific cases of Mk (fc = and 1, respectively). We consider the k-th order 
weighted mean fragment mass Mk+i/Mk, and assume the scaling, 



Mo 



(1) 



M — m min . , 
Ml 




(2) 



M fc+ i 



(3) 
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Figure 1. The maximum fragment mass M max vs. the imparted energy per unit sample mass e. 
Although the rough correlation between M max and e can be seen , it is too fluctuating to discuss 
quantitatively. 



In Fig. |3 we show the log- log plot of M2/M1 as a function of /i. Contrary to the 
Fig. ^ the data in Fig. [5] seem to be fitted by a unified scaling line. The scaling 
result for k = 1 is presented as a solid line in Fig. [21 We obtained the nontrivial 
scaling exponent a± — 0.84 ± 0.05. For other k regime, multi-scaling exponent 
values of <7fc were obtained as shown in Fig. [3] (circle marks). In spite of the trivial 
value <7fc = o = 1, ak varies with k, and seems to approach to the nontrivial value 
(-0.6). 



From the definition of 7^ as 



Mi M ' W 



the obtained jf. values are plotted as square marks in Fig. [3] It seems that 7^ 
approaches to the value around 0.6 again. Of course, Eqs. Q and (@J relate to 
each other. The relation Yli=i a i = holds for any k. Thus, when the er^ has 
a trivial value 1 for all fc, 7^ varies as (k — l)/fc. The trivial curves are shown as 
broken lines in Fig. [3] In addition, the relation (fc + 1)7^+1 — kjk — &k can be 
computed from Eqs. © and J3J|. If the difference between 7^+1 and jk becomes 
small, 7/c and <Jk approximately have the same value, as seen in Fig. for large k. 

On the other hand, when the imparted energy was sufficient to fully shatter, 
many fragments were created and the cumulative distribution of fragment mass 
obeyed power-law form. Our results suggest the power is 0.5 about the 2-D frag- 
mentation with the flat impact. In this regime, cumulative distribution functions 
are collapsed by the scaling function written as N(m)/Mo ~ f{m/(x~ rT ). The func- 
tion f(x) consists of the scaling part f(x) ~ x~ a5 and the decaying part due to the 
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Figure 2. The weighted mean fragment mass M2 /Mi vs. the multiplicity fi. The solid line indicates 
the form of the power-law fitting M2/M1 ~ n~ ai (a± = 0.84 ± 0.05). Three different size results 
are fitted by the unified scaling independently of size. 




Figure 3. Multi-scaling exponent <r fc and 7^ obtained from fc-th order weighted mean fragment 
mass. The broken lines indicate the case of trivial integral scaling corresponding to m = 1 and 
Ik = (k — f)/fc- The solid lines depict results of the binomial multiplicative model at a = 2/3. 
The cr fe and asymptotically approach to the same value, in large k range. From the definition 
of Eqs. an d 0i the values ct^—q and 7fe=l exactly show f and 0, respectively. 
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finite size effect. The value 0.5 concurs to the results of Hayakaw£p>2l and Astrom et 
alP. In contrast, this value is not consistent with the percolation scaling ansatzE] 
In the percolation scaling ansatz, the scaling exponent of cluster size cumulative 
distribution must be greater than 1. Therefore, we can consider that the univer- 
sality classes of critical fragmentation and percolation criticality are different each 
other. 



3 Model 

In order to explain this multi-scaling property, we introduce a simple biased cascade 
model. A binomial multiplicative process is considered with a unit mass initial 
condition. Here we consider a asymmetrical cleaving presented by a parameter a. 
We can limit the range of the parameter oasl/2<a<lby the symmetry of the 
model. The initial unit mass is divided into two fragments whose masses are a and 
1 — a at first step. This biased cleaving continues some steps until the imparted 
energy dissipates. In this model, we can easily calculate the exponents au and 7& 
from Eqs. © and as 

a k + (l-a) k v ; 

a k + (l-a) k = 2"^*, (5b) 

or more explicit forms as, 

ln[a fe+1 + (1 - a) k+1 ] - ln[a fe + (1 - a) k ] 
o k = — , (6a) 

If we choose a value a — 2/3, the a k and 7^ become the values depicted by the solid 
lines in Fig. EI O ne can confirm the pretty good agreement with experimental data. 
The trivial case presented by broken lines in Fig.[3]corresponds to the case a = 1/2. 
In this case, all fragments at each step are perfectly equal. In the case a 7^ 1/2, the 
fragment size distribution has variation and exhibits multifractal scaling. 

This model is so simple that we can calculate the fragment mass and the number 
of fragments exactly. We consider the s-th step, and suppose the fragments in which 
the factor a (or 1 — a) works t (or s — t) times. In such fragments, the mass m s (t) 
is written as, 

m s (t)=a t (l-a) s -\ (i < a < 1). (7) 
And the number of fragments n s (t) is described as, 



BinomialModel: submitted to World Scientific on February 2, 2008 5 




Figure 4. The cumulative fragment mass distribution for the binomial multiplicative model with 
a = 2/3 and s = 10. 



Since we are interested in the cumulative distribution of fragment mass, the 
cumulative number of fragments N s (t) is calculated as 



N s (t) 



/'CO 

/ n s (t')dt' = J2^ 



(9) 



In Fig. 01 we show the cumulative distribution of fragment mass obtained from 
the model. The parameters are taken as a = 2/3 and s = 10. The line of slope 
—0.5 corresponding to the experimental result is also shown as a solid line in Fig. 
01 Unfortunately, clear power-law, which follows the experimental data, cannot be 
observed. However, the distribution curve in Fig. 0] seems to include the region of 
slope —0.5. In this model, we can calculate the local power r — 1 directly by the 
relation, 



N a {t-1) 
N s (t) 



m s (t - 1) 
m s (t) 



(10) 



Solving the Eq. . we obtain the exact form of r — 1 as follows, 

Es s! 
i=t-l i\(s-i)\ 



T - 1 = -- 



In 


~N s (t-iy 

N s (t) 


In 
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m s (t - 1) 
m s (t) 


In 



E 



i=t i\(s—i)\ 



,t-l 



(1 



-t+1 



a*(l - a) 6 



(11) 



We show the relations among r — 1, s, and t in Fig. |SJ As can be seen in Fig. [31 
the lower limit of the local slope r — 1 is 0, and it has a divergent tendency. The 
value 0.5 is not a particular one. 
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Figure 5. The local scaling exponent of cumulative fragment mass distribution r — 1. (a) The 
relations between r — 1 and t are shown. Each curve corresponds to the case s = 5, 6, ■ • • , 10 from 
left to right, (b) The relations between r — 1 and s are shown. Each curve corresponds to the 
case t = 5, 6, • • • , 10 from left to right. 



4 Discussion 

In Sec. we concluded that the universality of the critical fragmentation differs 
from that of percolation. Instead, the weighted mean fragment mass was studied to 
reveal the universality of the critical fragmentation. It indicates the multi-scaling 
nature and is modeled by the simple binomial multiplicative model. There are 
some other candidates for the critical fragmentation. From now on, we discuss and 
compare them. 

Similar multiplicative model for turbulent flows was proposed by Meneveau and 
breemvasa rP. They investigated the energy cascade of eddies, and obtained good 
agreement with experimental data at a — 0.3 (in their paper, the corresponding 
parameter was written as p\). This value slightly coincides to ours 1 — a = 1/3. 
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The same physical mechanism might dominate both cascades of fragmentation and 
turbulence. 

We can fit the data by a = 2/3 very well indeed, however, the reason of symme- 
try breaking by a ^ 1/2 is not understood well. While the model always requires 
the exact asymmetry presented by a, the cleaving point might distribute. We can 
consider the simple distributed model as described below. We set the unit mass 
segment initial condition again, and consider the probability p(x)dx which presents 
the cleaving point in the range from x to x + dx at each step. We assume the 
symmetrical distribution as p(x) = Ax (0 < x < 1/2), and 4 — Ax (1/2 < x < 1). 
This is one of the simplest distribution presented by isosceles triangles. The nor- 
malization condition of this model is Jq p(x)dx = 1. In this model, we can calculate 
the expectation value of the cleaving point x (or 1 — x) as, 



r 



(1 - x)Axdx + x(A- Ax)dx = -. (12) 



Note that we cannot distinguish the cleaving state (x, 1 — x) and the state(l — x, x). 
Thus, the x can be limited in the region 1/2 < x < 1. The expectation value is 
nearly the same as one (a = 2/3) of the above mentioned multiplicative model. We 
can also calculate the fc-th order moment M k as, 



M k = / [x k + (1 - x) k ]p(x)dx = 



o 1 v ' J ' w (fc + l)(fc + 2) 



I \ fe+i' 



(13) 



We show the a k computed from the Eqs. I|13l) . Q, and P as a solid line in Fig. EI 
The result does not supply the agreement with the experimental data, particularly 
in large k range. 

Matsuhitcdland Turcottd"^ introduced the model for power-law fragmentation. 
Matsushita examined the model in which each fragment cleaves into 4 pieces at each 
step. Then 1 piece does not break any more, and the other 3 pieces cleave into 4 
sub-pieces at next step. The same procedure works upon all sub-pieces at each step. 
According to his model, the exponent of the power-law cumulative distribution of 
fragment mass becomes r — 1 = In3/ln4 ~ 0.79. We can easily modify this model 
to the case t — 1 = In 2/ In 4 = 1/2 by changing the remaining piece number 1 into 
2. In this condition, we can calculate the Mk for this modified version of the partial 
remaining model at s-th step as, 

s / , lk 



i=l v 7 

Then, the <y% can be computed again, however, the value of depends not only 
on k, but also on s. We show the a k obtained from this model at s = 10 as a 
broken line in Fig. [fj] This model also cannot provide the appropriate curve of the 
CTfc. Thus the exact a = 2/3 binomial multiplicative model seems to be the most 
possible model in terms of multi-scaling exponents a k . 
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Figure 6. Comparison between the experimental data and the other considerable models. The 
experimental data are shown as circle marks. The exponent ay. obtained from the symmetric 
(isosceles triangle) distribution model is presented by the solid line. The broken line indicates the 
result of the modified partial remaining model at s = 10. 



5 Conclusions 

We investigated the criticality of brittle fragmentation. Some models to interpret 
the experimental result are proposed. The exact binomial multiplicative model can 
produce the adequate approximation for the exponent <Jk ■ And the cumulative dis- 
tribution obtained from the model is not so worse. However, it requires that the 
cleaving point is exactly at a. Since the isosceles triangle model has non-divergent 
standard deviation, the distribution of fragment mass resulting from the m odel m ust 
approach to the log-normal distribution due to the central limit theoreti ca l In 
addition, its Ok differs from the experimental data, particularly in large k region. 
The modified partial remaining model can explain the experimental value of the ex- 
ponent t very well. However, the <ik from the model shows large discrepancy from 
the experimental data. Each model has merits and demerits. The totally sufficient 
model is not presented yet. Furthermore, these scaling exponents will depend on 
detail load condition and dimensionality of fractured object. More detailed exper- 
iments and analyses are necessary to solve the criticality of brittle fragmentations 
completely. 
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